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Abstract 



A fundamental goal in QCD is to understand the non-perturbative structure of 
hadrons at the amplitude level — not just the single-particle flavor, momentum, and 
helicity distributions of the quark constituents, but also the multi-quark, gluonic, and 
hidden-color correlations intrinsic to hadronic and nuclear wavef unctions. A natural 
calculus for describing the bound-state structure of relativistic composite systems in 
quantum field theory is the light-front Fock expansion which encodes the properties of 
a hadrons in terms of a set of frame-independent n— particle wavefunctions. Light-front 
quantization in the doubly-transverse light-cone gauge has a number of remarkable ad- 
vantages, including explicit unitarity, a physical Fock expansion, the absence of ghost 
degrees of freedom, and the decoupling properties needed to prove factorization theo- 
rems in high momentum transfer inclusive and exclusive reactions. A number of ap- 
plications are discussed in these lectures, including semileptonic B decays, two-photon 
exclusive reactions, and deeply virtual Compton scattering. The relation of the intrinsic 
sea to the light-front wavefunctions is discussed. A new type of jet production reaction, 
"self-resolving diffractive interactions" can provide direct information on the light-front 
wavefunctions of hadrons in terms of their quark and gluon degrees of freedom as well 
as the composition of nuclei in terms of their nucleon and mesonic degrees of freedom. 



1 Introduction 

In principle, quantum chromodynamics provides a fundamental description of hadron 
and nuclear physics in terms of quark and gluon degrees of freedom. QCD has been 
developed and successfully tested extensively, particularly in inclusive and exclusive 
processes involving collisions at large momentum transfer where factorization theorems 
and the smalhicss of the QCD effective coupling allow perturbativc predictions. How- 
ever, despite its empirical successes, many fundamental questions about QCD have not 
been resolved. These include a rigorous proof of color confinement, the behavior of the 
QCD coupling at small momenta, a rigorous description of the structure of hadrons 
in terms of their quark and gluon degrees of freedom, the problem of asymptotic n! 
growth of the perturbation theory (rcnormalon phenomena), the nature of the pomeron 
and Reggeons, the nature of shadowing and ant i- shadowing in nuclear collisions, the 
apparent conflict between QCD vacuum structure and the small size of the cosmolog- 
ical constant, and the problems of scale and scheme ambiguities in perturbative QCD 
expansion. 

In these lectures I will focus on one of the central questions in QCD — the non- 
perturbative description of the proton and other hadrons as composites of confined. 
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relativistic quark and gluon quanta. The goal is a frame-independent, quantum- 
mechanical representation of hadrons at the amplitude level capable of encoding multi- 
quark, hidden-color and gluon momentum, helicity, and flavor correlations in the form 
of universal process-independent hadron wavefunctions. For example, the measurement 
and interpretation of the basic parameters of the electroweak theory and CP violation 
depends on an understanding of the dynamics and phase structure of B decays at the 
amplitude level. As I will discuss in these lectures, light-front quantization allows a 
unifying wavefunction representation of non-perturbative hadron dynamics in QCD. 
Remarkably, it is now possible to measure the wavefunctions of a relativistic hadron by 
diffractively dissociating it into jets whose momentum distribution is correlated with 
the valence quarks' momenta |^, ^, 3. It is also particularly important to under- 
stand the shape of the gauge- and process-independent meson and baryon valence-quark 
distribution amplitudes 4>m{x,Q), and (f)B{xi,Q). These quantities specify how a 
hadron shares its longitudinal momentum among its valence quarks; they control vir- 
tually all exclusive processes involving a hard scale Q, including form factors, Compton 
scattering and photoproduction at large momentum transfer, as well as the decay of a 
heavy hadron into specific final states P, 0. 

What do we mean by a hadronic wavefunction, and what are its degrees of freedom? 
In the chiral effective Lagrangian approximation, a baryon can be represented as a 
classical soliton solution as in the Skyrme model. On the other hand, bag models, the 
observed baryon spectroscopy, and magnetic moment phenomenology suggests that 
baryons are composites of three "constituent" quarks or perhaps a quark bound to an 
effective spin-0 or spin-1 diquark. The cloudy bag model takes into account the 
effect of meson-baryon fluctuations such as mr'^ and AK~^ in the proton. This in turn 
implies the existence of sea quarks with an u{x) 7^ d{x) and s{x) 7^ s{x) asymmetries 
in the momentum and spin distributions of the quark sea. 

The ra-particle Schrodinger momentum space wavefunction ipN{Pi) of a composite 
system is the projection of the exact eigenstate of the equal-time Hamiltonian on the 
n-particle states of the non-interacting Hamiltonian, the Fock basis. It represents the 
amplitude for finding the constituents with three-momentum pi, orbital angular mo- 
mentum, and spin, subject to three-momentum conservation and angular momentum 
sum rules. The constituents are on their mass shell, Ei = ^Jpf + m? but do not con- 
serve energy ^2^=1 > E = y/p^ + M'^. However, in a relativistic quantum theory, 
a bound-state cannot be represented as a state with a fixed number of constituents. 
For example, the existence of gluons which propagate between the valence quarks nec- 
essarily implies that the hadron wavefunction must describe states with an arbitrary 
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number of gluons. Thus a hadronic wavefunction must describe fluctuations in particle 
number n, as well as momenta and spin. One has to take into account fluctuations in 
the wavefunction which allow for any number of sea quarks, as long as the total quan- 
tum numbers of the constituents are compatible with the overall quantum numbers of 
the baryon. 

Gluon exchange between quarks of different nucleons in a nucleus will leave the 
nucleons in an exited octet state of SU{3)c, thus requiring the existence of "hidden 
color" configurations in a nuclear wavefunction. What is worse is that the description 
of the state at a given time to depends not only on the choice of gauge but also on 
the Lorentz frame of the observer since two observers will differ on their definition of 
simultaneity. 

A number of non-perturbative formalisms have been developed which in principle 
could provide the exact spectra and wavefunctions of bound states in relativistic quan- 
tum field theory, including QCD. These include Euclidean lattice gauge theory, the 
covariant Bethe-Salpeter bound state equation, and Hamiltonian methods based on 
equal-time or light-front quantization. Lattice gauge theory, particularly in the heavy 
quark approximation, has been an enormously helpful guide to the hadronic and glue- 
ball spectrum of QCD; however, progress in evaluating hadron wavefunctions has been 
so far limited to the evaluation of certain moments of hadronic matrix elements f^. 

The Bethe-Salpeter formalism in principle provides a rigorous and systematic, renor- 
malizable, and explicitly covariant treatment of bound states of two interacting fields. 
There has been recent progress applying model forms of Dyson- Schwinger equations to 
the Bethe-Salpeter formalism in order to incorporate the effects of dynamical symme- 



try breaking in QCD [10]. However, the actual interaction kernel in the bound state 
equation and the evaluation of current matrix elements requires a sum over an infinite 
number of irreducible graphs; the use of a single kernel (the ladder approximation) 
breaks gauge invariance, crossing symmetry, and has an improper heavy quark limit. 
It is also intrinsically difficult to apply the Bethe Salpeter to bound states of three or 
more constituents. 

There has been progress applying equal-time Hamiltonian methods to QCD in 



Coulomb gauge |TT[; however, the renormalization procedure in this formalism is not 
well-defined and the results are typically frame- dependent. The evaluation of current 
matrix elements in the equal-time formalism is especially problematic since one must 
include contributions to the current from vacuum excitations. For similar reasons, the 
expansion of an eigenstate on the equal-time Fock basis of free quanta is not well- 
defined. 



4 



In the light-front quantization, one takes the hght-cone time variable t + z/c as 
the evolution parameter instead of ordinary time t. (The z direction is an arbitrary 
reference direction.) The method is often called "light-front" quantization rather than 
"light-cone" quantization since the equation a;"^ = r = defines a hyperplane corre- 
sponding to a light-front. The light-front fixes the initial boundary conditions of a 
composite system as its constituents are intercepted by a light-wave evaluated at a 
specific value oi = t + z / c. In contrast, determining an atomic wavef unction at a 
given instant t = to requires measuring the simultaneous scattering of Z photons on 
the Z electrons. 

An extensive review and guide to the light-front quantization literature can be found 
in Ref. [|12|]. I will use here the notation = A—, A±), where = A^ ±A^ and 



the metric A- B = UA+B- + A-B+) - A^ ■ B^. 



The origins of the light-front quantization method can be traced to Dirac |]T3| who 
noted that 7 out of the 10 Poincare' generators, including a Lorentz boost K^, are 
kinematical (interaction-independent) when one quantizes a theory at fixed light-cone 
time. This in turn leads to the remarkable property that the light-front wavefunctions 
of a hadron are independent of the hadron's total momentum, whether it is at rest or 
moving. Thus once one has solved for the light-front wavefunctions, one can compute 
hadron matrix elements of currents between hadronic states of arbitrary momentum. 
In contrast, knowing the rest frame wavefunction at equal time, does not determine 
the moving hadron's wavefunction. 

Light-front wavefunctions are related to momentum-space Bethe-Salpeter wavefunc- 
tions by integrating over the relative momenta k~ = — since this projects out 

= 0. The light-front Fock space is the eigenstates of the free light-front Hamiltonian; 
i.e., it is a Hilbert space of non-interacting quarks and gluons, each of which satisfy 
fc^ = and k~ = ^ > 0. The light-front wavefunctions are the projections of the 
hadronic eigenstate on the light-front Fock basis, the complete set of color singlet states 
of the free Hamiltonian. An essentially equivalent approach, pioneered by Weinberg 
T\ IT5[, is to evaluate the equal-time theory from the perspective of an observer moving 



in the negative z direction with arbitrarily large momentum — * — oo. The light-cone 
fraction x = ^ of a constituent can be identified with the longitudinal momentum 
X — in a hadron moving with large momentum P^. 

It is convenient to define the invariant light-front Hamiltonian: H^c^ = P^ P^ ~ P]_ 
where P"^ = P^ ± P^. The operator P~ = ij;: generates light-cone time transla- 
tions. The P~^ and Pi momentum operators are independent of the interactions, and 
thus are conserved at all orders. The eigen-spectrum of H^c^ in principle gives the 
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entire mass squared spectrum of color-singlet hadron states in QCD, together with 
their respective light-front wavef unctions. For example, the proton state satisfies: 
H^Q^l \E'p) = Mp \ \E'p). The projection of the proton's eigensolution | \E'p) on the color- 
singlet B = 1, Q = 1 eigenstates {| ^)} of the free Hamiltonian H^c^ig = 0) gives the 
light-front Fock expansion: Wi 



n>3,Ai •' ^/x'ilQ^T 



n 



fl] XiP~^ , XiP± + k±i, Xi^ k±i, Aj). 

It is especially convenient to develop the light-front formalism in the light-cone 
gauge A'^ = + = 0. In this gauge the A~ field becomes a dependent degree of 
freedom, and it can be eliminated from the gauge theory Hamiltonian, with the addition 
of a set of specific instantaneous light-cone time interactions. In fact in QCD{1 + 
1) theory, this instantaneous interaction provides the confining linear x~ interaction 
between quarks. In 3 + 1 dimensions, the transverse field A-^ propagates massless spin- 
one gluon quanta with polarization vectors [Q which satisfy both the gauge condition 

= and the Lorentz condition k ■ e = 0. Thus no extra condition on the Hilbert 
space is required. 

The light-front Fock wavefunctions ipn/nixi, k±i, Xi) thus interpolate between the 
hadron H and its quark and gluon degrees of freedom. The light-cone momentum frac- 
tions of the constituents, = kf /P^ with X^^Li^i — 1; the transverse momenta 
k±i with = 0_L appear as the momentum coordinates of the light-front Fock 

wavefunctions. A crucial feature is the frame-independence of the light-front wavefunc- 
tions. The Xi and k±i are relative coordinates independent of the hadron's momentum 
P^. The actual physical transverse momenta are = XiP± + k±i. The A, label the 
light-front spin 5*^ projections of the quarks and gluons along the z direction. The 
physical gluon polarization vectors e^{k, X = ±1) are specified in light-cone gauge by 
the conditions k ■ e = 0, rj ■ e = e'^ = 0. Each light-front Fock wavefunction satis- 
fies conservation of the z projection of angular momentum: = Yh=i + S"=i ■ 
The sum over represents the contribution of the intrinsic spins of the n Fock state 
constituents. The sum over orbital angular momenta /| = —i{kj^ — k'j-^) derives 
from the n — 1 relative momenta. This excludes the contribution to the orbital angular 
momentum due to the motion of the center of mass, which is not an intrinsic property 



of the hadron ITT 
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The interaction Hamiltonian of QCD in light-cone gauge can be derived by system- 



atically applying the Dirac bracket method to identify the independent fields [|T8|. It 
contains the usual Dirac interactions between the quarks and gluons, the three-point 
and four-point gluon non-Abelian interactions plus instantaneous light-front-time gluon 
exchange and quark exchange contributions 



-.1 



where 



I 3 jabc jade A'^^A A^'^ 

-—7+ ^ ?+ (1) 
2 " (9 )^ " 



J^a = i^Vita^^' + faUd-AM"' ■ (2) 



In light-cone time-ordered perturbation theory, a Green's functions is expanded as 
a power series in the interactions with light-front energy denominators I^initiai k~ — 
^intermediate + replacing the usual energy denominators. [For a review see Ref. 



T6| .] In general each Feynman diagram with n vertices corresponds to the sum of n\ 
time-ordered contributions. However, in light-cone-time-ordered perturbation theory, 
only those few graphs where all kf > survive. In addition the form of the light-front 



kinetic energies is rational: k~ = , replacing the nonanalytic = yk"^ + m? of 

equal-time theory. Thus light-cone-time-ordered perturbation theory provides a viable 
computational method where one can trace the physical evolution of a process. The 
integration measures are only three-dimensional d'^k±dx; in effect, the k~ integral of 
the covariant perturbation theory is performed automatically. 

Alternatively, one derive Feynman rules for QCD in light-cone gauge, thus allow- 
ing the use of standard covariant computational tools and renormalization methods 
including dimensional regularization. Prem Srivastava and I have recently pre- 



sented a systematic study of light-front-quantized gauge theory in light-cone gauge 
using a Dyson- Wick S-matrix expansion based on light- cone-time-ordered products. 
The gluon propagator has the form 

0| TiA%ix)A\m |0) = 1^ / d^k e--- §^ (3) 
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where we have defined 



Here ra^ is a null four-vector, gauge direction, whose components are chosen to be 
= 6^~^, = 5'^-. Note also 

D^x{k)D\{k) = D^^{k)D^,{k) = -D^,{k), 

k^D^,{k) = 0, n^D^,{k) = D^,{k) = 0, 

D,^{q) D^^{k) D^^iq') = -Dx^{q)D^P{q'). (5) 

The gauge field propagator i D ^^^ik) / iJxP' + ie) is transverse not only to the gauge 
direction but also to fc^, i.e., it is doubly-transverse. This leads to appreciable 
simplifications in the computations in QCD. For example, the coupling of gluons to 
propagators carrying high momenta is automatic. The absence of coUinear divergences 
in irreducible diagrams in the light-cone gauge greatly simplifies the leading-twist fac- 
torization of soft and hard gluonic corrections in high momentum transfer inclusive 
and exclusive reactions since the numerators associated with the gluon coupling 
only have transverse components. The renormalization factors in the light-cone gauge 
are independent of the reference direction n^. Since the gluon only has physical polar- 
ization, its renormalization factors satisfy Zi = Z3. Because of its explicit unitarity in 
each graph, the doubly-transverse gauge is well suited for calculations identifying the 
"pinch" effective charge |T9|, pO |. 



The running coupling constant and QCD function have also been computed at 
one loop in the doubly-transverse light-cone gauge |]TB|. It is also possible to effectively 
quantize QCD using light- front methods in covariant Feynman gauge 

A remarkable advantage of light-front quantization is that the vacuum state | 0) of 
the full QCD Hamiltonian evidently coincides with the free vacuum. The light-front 
vacuum is effectively trivial if the interaction Hamiltonian applied to the perturbative 
vacuum is zero. Note that all particles in the Hilbert space have positive energy = 
^{k~^ + k~), and thus positive light-front fc^. Since the plus momenta J2 kf is conserved 
by the interactions, the perturbative vacuum can only couple to states with particles in 
which all kf = 0; i.e., so called zero-mode states. In the case of QED, a massive electron 
cannot have k~^ = unless it also has infinite energy. In a remarkable calculation. 



Bassetto and collaborators [|22[ have shown that the computation of the spectrum of 
QCD[1 -|- 1) in equal time quantization requires constructing the full spectrum of non 
perturbative contributions (instantons). In contrast, in the light-front quantization of 
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gauge theory, where the = singularity of the instantaneous interaction is defined 
by a simple infrared regularization, one obtains the correct spectrum of QCD{1 + 1) 
without any need for vacuum-related contributions. 

In the case of QCD(3+1), the momentum-independent four-gluon non-Abelian in- 
teraction in principle can couple the perturbative vacuum to a state with four collinear 
gluons in which all of the gluons have all components /cf = 0, thus hinting at role for 
zero modes in theories with massless quanta. In fact, zero modes of auxiliary fields 
are necessary to distinguish the theta- vacua of massless QED(1+1) |2^, or to 

represent a theory in the presence of static external boundary conditions or other con- 
straints. Zero-modes provide the light-front representation of spontaneous symmetry 
breaking in scalar theories ||2^ . 



There are a number of other simplifications of the light-front formalism: 

1. The light-front wavefunctions describe quanta which have positive energy, pos- 
itive norm, and physical polarization. The formalism is thus physical, and unitary. 
No ghosts fields appear explicitly, even in non-Abelian theory. The wavefunctions are 
only functions of three rather than four physical momentum variables: the light-front 
momentum fractions Xi and transverse momenta fc^. The quarks and gluons each have 
two physical polarization states. 

2. The set of light-front wavefunctions provide a frame-independent, quantum- 
mechanical description of hadrons at the amplitude level capable of encoding multi- 
quark and gluon momentum, helicity, and fiavor correlations in the form of universal 
process-independent hadron wavefunctions. Matrix elements of spacelike currents such 
as the spacelike electromagnetic form factors have an exact representation in terms of 
simple overlaps of the light-front wavefunctions in momentum space with the same Xi 
and unchanged parton number |2^, In the case of timelike decays, such as 
those determined by semileptonic B decay, one needs to include contributions in which 
the parton number An = 2. 0]. The leading- twist off- forward parton distributions 
measured in deeply virtual Compton scattering have a similar light-front wavefunction 
representation |3^. 

3. The high x ^ 1 and high limits of the hadron wavefunctions control pro- 
cesses and reactions in which the hadron wavefunctions are highly stressed. Such con- 
figurations involve far-off-shell intermediate states and can be systematically treated in 
perturbation theory |3^, |^ . 

4. The leading-twist structure functions qi{x,Q) and g{x,Q) measured in deep 
inelastic scattering can be computed from the absolute squares of the light-front wave- 
functions, integrated over the transverse momentum up to the resolution scale Q. All 
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helicity distributions are thus encoded in terms of the hght-front wavefunctions. The 
DGLAP evolution of the structure functions can be derived from the high k±_ prop- 
erties of the hght-front wavefunctions. Thus given the hght-front wavefunctions, one 
can compute aU of the leading twist helicity and transversity distributions measured 
in polarized deep inelastic lepton scattering. For example, the helicity-specific quark 
distributions at resolution A correspond to 

gVA.(a;,A) = E / n^^TS^El<k^.,^^.,A.)p (6) 




Six - x,)S,,xM^' - ^l) ' 

where the sum is over all quarks qa which match the quantum numbers, light-front 
momentum fraction x, and helicity of the struck quark. Similarly, the transversity 
distributions and off-diagonal helicity convolutions are defined as a density matrix of 
the light-front wavefunctions. This defines the LC factorization scheme where the 
invariant mass squared = J27=i {ki.i + of the n partons of the light-front 

wavefunctions is limited to < A^ 

5. The distribution of spectator particles in the final state in the proton fragmen- 
tation region in deep inelastic scattering at an electron-proton collider are encoded in 
the light-front wavefunctions of the target proton. Conversely, the light-front wave- 
functions can be used to describe the coalescence of comoving quarks into final state 
hadrons. 

6. The light-front wavefunctions also specify the multi-quark and gluon correlations 
of the hadron. Despite the many sources of power-law corrections to the deep inelastic 
cross section, certain types of dynamical contributions will stand out at large Xbj since 
they arise from compact, highly-correlated fluctuations of the proton wavefunction. 
In particular, there are particularly interesting dynamical (9(1/(5^) corrections which 
are due to the interference of quark currents; i.e., contributions which involve leptons 



scattering amplitudes from two different quarks of the target nucleon |34 . 

7. The higher Fock states of the light hadrons describe the sea quark structure 
of the deep inelastic structure functions, including "intrinsic" strangeness and charm 
fluctuations specific to the hadron's structure rather than gluon substructure ||35|, |36[| . 
Ladder relations connecting state of different particle number follow from the QCD 
equation of motion and lead to Regge behavior of the quark and gluon distributions at 
X |7|. 

8. The gauge- and process-independent meson and baryon valence-quark distribu- 
tion amplitudes (f)M{x,Q), and 4>B{xi,Q) which control exclusive processes involving a 
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hard scale Q, including heavy quark decays, are given by the valence light-front Fock 
state wavefunctions integrated over the transverse momentum up to the resolution scale 
Q. The evolution equations for distribution amplitudes follow from the perturbative 
high transverse momentum behavior of the light-front wavefunctions |[T6| . 

9. The line-integrals needed to defining distribution amplitudes and structure func- 
tions as gauge invariant matrix elements of operator products vanish in light-front 
gauge. 

10. Proofs of factorization theorems in hard exclusive and inclusive reactions are 
greatly simplified since the propagating gluons in light-cone gauge couple only to trans- 
verse currents; coUinear divergences are thus automatically suppressed. 

11. At high energies each light-front Fock state interacts distinctly; e.g., Fock 
states with small particle number and small impact separation have small color dipole 
moments and can traverse a nucleus with minimal interactions. This is the basis for 
the predictions for "color transparency" in hard quasi-exclusive [p8t |39| and diffractive 
reactions [|, |, |] . 

12. The Fock state wavefunctions of hadron can be resolved by a high energy 
diffractive interaction, producing forward jets with momenta which follow the light- 
front momenta of the wavefunction 0, ^ 

13. The deuteron form factor at high is sensitive to wavefunction configura- 
tions where all six quarks overlap within an impact separation b±i < 0{1/Q). The 
leading power-law fall off predicted by QCD is F^iQ'^) = /(as(<5^))/(Q^)^, where, 
asymptotically, /(as(<5^)) oc as(Q^)^^^'^ 0! EH- general, the six-quark wavefunc- 
tion of a deuteron is a mixture of five different color-singlet states. The dominant color 
configuration at large distances corresponds to the usual proton-neutron bound state. 
However at small impact space separation, all five Fock color-singlet components even- 
tually evolve to a state with equal weight, i.e., the deuteron wavefunction evolves to 
80% "hidden color" |^T|. The relatively large normalization of the deuteron form fac- 
tor observed at large hints at sizable hidden-color contributions Hidden color 
components can also play a predominant role in the reaction 'yd J / ippn at threshold 
if it is dominated by the multi-fusion process 'jgg J/ip Hard exclusive nuclear 
processes can also be analyzed in terms of "reduced amplitudes" which remove the 
effects of nucleon substructure. 

Light-front wavefunctions are thus the frame-independent interpolating functions 
between hadron and quark and gluon degrees of freedom. Hadron amplitudes are 
computed from the convolution of the light-front wavefunctions with irreducible quark- 
gluon amplitudes. More generally, all multi-quark and gluon correlations in the bound 
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state are represented by the light-front wavefunctions. The hght-front Fock represen- 
tation is thus a representation of the underlying quantum field theory. I will discuss 
progress in computing light-front wavefunctions directly from QCD in Sections 9 and 
10. 



2 Other Theoretical Tools 

In addition to the light-front Fock expansion, a number of other useful theoretical tools 
are available to eliminate theoretical ambiguities in QCD predictions: 

(1) Conformal symmetry provides a template for QCD predictions leading to 



relations between observables which are present even in a theory which is not scale 
invariant. For example, the natural representation of distribution amplitudes is in 
terms of an expansion of orthonormal conformal functions multiplied by anomalous 



dimensions determined by QCD evolution equations Thus an important 

guide in QCD analyses is to identify the underlying conformal relations of QCD which 
are manifest if we drop quark masses and effects due to the running of the QCD 
couplings. In fact, if QCD has an infrared fixed point (vanishing of the Cell Mann-Low 
function at low momenta), the theory will closely resemble a scale-free conformally 
symmetric theory in many applications. 



(2) Commensurate scale relations ^ are perturbative QCD predictions which 



relate observable to observable at fixed relative scale, such as the "generalized Crewther 



relation" which connects the Bjorken and Gross-Llewellyn Smith deep inelastic 
scattering sum rules to measurements of the e~^e~ annihilation cross section. Such 
relations have no renormalization scale or scheme ambiguity. The coefficients in the 
perturbative series for commensurate scale relations are identical to those of conformal 
QCD; thus no infrared renormalons are present One can identify the required 



conformal coefficients at any finite order by expanding the coefficients of the usual 
PQCD expansion around a formal infrared fixed point, as in the Banks-Zak method ||20|| . 
All non-conformal effects are absorbed by fixing the ratio of the respective momentum 
transfer and energy scales. In the case of fixed-point theories, commensurate scale 
relations relate both the ratio of couplings and the ratio of scales as the fixed point is 
approached . 



(3) ay and Skeleton Schemes. A physically natural scheme for defining the QCD 
coupling in exclusive and other processes is the ay(Q^) scheme defined from the poten- 
tial of static heavy quarks. Heavy-quark lattice gauge theory can provide highly precise 
values for the coupling. All vacuum polarization corrections due to fermion pairs are 
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then automatically and analytically incorporated into the Gell Mann-Low function, 
thus avoiding the problem of explicitly computing and resumming quark mass correc- 
tions related to the running of the coupling The use of a finite effective charge 
such as ay as the expansion parameter also provides a basis for regulating the infrared 
nonperturbative domain of the QCD coupling. A similar coupling and scheme can be 
based on an assumed skeleton expansion of the theory pO[] . 

(4) The Abelian Correspondence Principle. One can consider QCD predictions as 
analytic functions of the number of colors and flavors Np. In particular, one can 
show at all orders of perturbation theory that PQCD predictions reduce to those of an 
Abelian theory at A'',^ — with S = CpOi^ and Np = 2Np / Cp held fixed |Q . There is 
thus a deep connection between QCD processes and their corresponding QED analogs. 



3 Applications of light-front Wavefunctions to Cur- 
rent Matrix Elements 

The light-front Fock representation of current matrix elements has a number of simpli- 
fying properties. The space-like local operators for the coupling of photons, gravitons 
and the deep inelastic structure functions can all be expressed as overlaps of light-front 
wavefunctions with the same number of Fock constituents. This is possible since one 
can choose the special frame g"*" = ||2^, ^ for space-like momentum transfer and take 
matrix elements of "plus" components of currents such as and No contribu- 

tions to the current matrix elements from vacuum fluctuations occur. Similarly, given 
the local operators for the energy-momentum tensor T^^{x) and the angular momentum 
tensor M'^'^^(x), one can directly compute momentum fractions, spin properties, and 
the form factors A{q^) and B{q^) appearing in the coupling of gravitons to composite 



systems [17 



In the case of a spin-i composite system, the Dirac and Pauli form factors 
and F2{q^) are defined by 



{P'\J^mP) = u{P') F^iq')r + F,iq 



2M 



(7) 



where q^ = (P' — P)^ and u{P) is the bound state spinor. In the light-front formalism 
it is convenient to identify the Dirac and Pauli form factors from the helicity-conserving 
and helicity-fiip vector current matrix elements of the current p9[] : 



J+(0) 



2P^ 




Fiiq' 



(8) 
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P + qA 



J+(0) 



2P^ 



2M 



(9) 



The magnetic moment of a composite system is one of its most basic properties. The 
magnetic moment is defined at the — > hmit, 



2M 



[Fi(0)+F2(0)] 



(10) 



where e is the charge and M is the mass of the composite system. We use the standard 
hght-front frame (g^ = g° ± q^): 



q = iq^,q ,q±) = 







-,q± 



M2 



where = —2P ■ q = — is 4-momentum square transferred by the photon. 

The Pauh form factor and the anomalous magnetic moment k, = 21? -^2(0) can then 
be calculated from the expression 



- {q^ - 



2, P2{q') 
2M 



'12) 



where the summation is over all contributing Fock states a and struck constituent 
charges ej. The arguments of the final-state light-front wavefunction are 



for the struck constituent and 



k'±i = k^i + (1 - Xi)qj 



k±i — k±i Xiq± 



(13) 



(14) 



for each spectator. Notice that the magnetic moment must be calculated from the spin- 
flip non-forward matrix element of the current. In the ultra-relativistic limit where the 
radius of the system is small compared to its Compton scale 1/M, the anomalous 
magnetic moment must vanish pOf. The light-front formalism is consistent with this 
theorem. 

The form factors of the energy- momentum tensor for a spin-^ composite are defined 

by 



(P'IT^'^(O)IP) = u{P') [A{q')j^^P'^ + B{q')—P''^a-'>-q, 

+C{q')^{q^q''-g^'''q')]uiP) 



(15) 



14 



where q^' = {P' - Py, P'' = \{P' + P)^, aH"^ = \{a^'h'' + a%^'), and u{P) is the 
spinor of the system. One can also readily obtain the light-front representation of the 
A{q^) and B{q^) form factors of the energy-tensor Eq. (|T3p |T7]. In the interaction 
picture, only the non- interacting parts of the energy momentum tensor T^'^{Q) need 
to be computed: 

T++(0) 



The and5(g2 
form factors Eqs. ( 
fraction x = k^/P^ 




2M 



(16) 
(17) 



2(P^ 

form factors Eqs. ([T6|) and (|T7|) are similar to the Fi{q^) and -^2(9^) 
]) and (P) with an additional factor of the light-front momentum 
of the struck constituent in the integrand. The B{q^) form factor 
is obtained from the non-forward spin-flip amplitude. The value of -B(O) is obtained in 
the — > limit. The angular momentum projection of a state is given by 



r 



X 



A(0) 



U) + [A{{}) + B{{))]u{P)-a'u{P) . (1^ 



This result is derived using a wave-packet description of the state. The {U) term is the 
orbital angular momentum of the center of mass motion with respect to an arbitrary 
origin and can be dropped. The coefficient of the {U) term must be 1; A(0) = 1 also 
follows when we evaluate the four- momentum expectation value (P^). Thus the total 
intrinsic angular momentum of a nucleon can be identified with the values of the 
form factors A{q^) and B{q'^) at = : 



1 . 

2 



[A(0) + B 



(19) 



The anomalous moment coupling P(0) to a graviton is shown to vanish for any 
composite system. This remarkable result, first derived by Okun and Kobzarev [^, 
55], ^ 0, is shown to follow directly from the Lorentz boost properties of the 



light-front Fock representation [17|. 

Dae Sung Hwang, Bo-Qiang Ma, Ivan Schmidt, and I [T?! have recently shown that 
the light-front wavefunctions generated by the radiative corrections to the electron in 
QED provides a simple system for understanding the spin and angular momentum 
decomposition of relativistic systems. This perturbative model also illustrates the in- 
terconnections between Fock states of different number. The model is patterned after 
the quantum structure which occurs in the one-loop Schwinger a/2TT correction to the 
electron magnetic moment In effect, we can represent a spin-| system as a com- 
posite of a spin- 1 fermion and spin-one vector boson with arbitrary masses. A similar 
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model has been used to illustrate the matrix elements and evolution of light-front helic- 



ity and orbital angular momentum operators |Q. This representation of a composite 
system is particularly useful because it is based on two constituents but yet is totally 
relativistic. We can then explicitly compute the form factors Fi{q'^) and -^2(0'^) of the 
electromagnetic current, and the various contributions to the form factors v4(g^) and 
B[q^) of the energy- momentum tensor. 
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Figure 1: Exact representation of electroweak decays and time-like form factors in the 
light-front Fock representation. 

Another remarkable advantage of the light-front formalism is that exclusive semilep- 
tonic -B-decay amplitudes such as -B — > Atv can also be evaluated exactly The 
time-like decay matrix elements require the computation of the diagonal matrix element 
n — > n where parton number is conserved, and the off-diagonal n+1 n—1 convolution 
where the current operator annihilates a qq' pair in the initial B wavefunction. See Fig. 
^ This term is a consequence of the fact that the time-like decay = {pi + > 
requires a positive light-front momentum fraction g"^ > 0. Conversely for space-like cur- 
rents, one can choose q~^ = 0, as in the Drell-Yan-West representation of the space-like 
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electromagnetic form factors. However, as can be seen from the explicit analysis of the 
form factor in a perturbative model, the off-diagonal convolution can yield a nonzero 
/q^ limiting form as g+ ^ 0. This extra term appears specifically in the case of 
"bad" currents such as J~ in which the coupling to qq fluctuations in the light-front 
wavefunctions are favored. In effect, the ^ limit generates 5{x) contributions as 
residues of the n + 1 ^ n — 1 contributions. The necessity for such "zero mode" 5{x) 
terms has been noted by Chang, Root and Yan [^, Burkardt [^, and Ji and Choi 

The off-diagonal n + 1 — > n — 1 contributions give a new perspective for the physics 
of 5-decays. A semileptonic decay involves not only matrix elements where a quark 
changes flavor, but also a contribution where the leptonic pair is created from the anni- 
hilation of a qq' pair within the Fock states of the initial B wavefunction. The semilep- 
tonic decay thus can occur from the annihilation of a nonvalence quark-antiquark pair 
in the initial hadron. This feature will carry over to exclusive hadronic -B-decays, such 
as —>■ iT~D~^. In this case the pion can be produced from the coalescence of a du 
pair emerging from the initial higher particle number Fock wavefunction of the B. The 
D meson is then formed from the remaining quarks after the internal exchange of aW 
boson. 

In principle, a precise evaluation of the hadronic matrix elements needed for B- 
decays and other exclusive electroweak decay amplitudes requires knowledge of all of 
the light-front Fock wavefunctions of the initial and final state hadrons. In the case of 
model gauge theories such as QCD(1+1) or coUinear QCD |^ in one-space and 
one-time dimensions, the complete evaluation of the light-front wavefunction is possible 
for each baryon or meson bound-state using the DLCQ method. It would be interesting 
to use such solutions as a model for physical S-decays. 



4 Light-front Representation of Deeply Virtual Comp- 
ton Scattering 

The virtual Compton scattering process ^(7*p — ^ 7p) for large initial photon virtu- 
ality = —Q^ has extraordinary sensitivity to fundamental features of the proton's 
structure. Even though the final state photon is on-shell, the deeply virtual process 
probes the elementary quark structure of the proton near the light cone as an effec- 
tive local current. In contrast to deep inelastic scattering, which measures only the 
absorptive part of the forward virtual Compton amplitude ImT^*p^^*p, deeply vir- 
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tual Compton scattering allows the measurement of the phase and spin structure of 
proton matrix elements for general momentum transfer squared t. In addition, the 
interference of the virtual Compton amplitude and Bethe-Heitler wide angle scattering 
Bremsstrahlung amplitude where the photon is emitted from the lepton line leads to 
an electron-positron asymmetry in the e^p —>■ e^jp cross section which is proportional 
to the real part of the Compton amplitude |6^, ^ The deeply virtual Compton 



amplitude 'y*p ^p is related by crossing to another important process 7*7 hadron 
pairs at fixed invariant mass which can be measured in electron-photon collisions ||67|| . 

To leading order in l/Q, the deeply virtual Compton scattering amplitude '~i*{q)p{P) 
'j{q')p{P') factorizes as the convolution in x of the amplitude t'^'^ for hard Compton scat- 
tering on a quark line with the generalized Compton form factors H{x, t, (), E{x, t, (), 
H{x, t, C), and E{x, t, Q of the target proton 0, |5|, |8|, |7g, |7l], 0, |7|, |7|, |75|, 0, [T^] . 



Here x is the light-front momentum fraction of the struck quark, and ( = Q^/2P ■ q 
plays the role of the Bjorken variable. The square of the four-momentum transfer from 
the proton is given by t = A = 2P ■ A = — -^^ — , where A is the difference of 
initial and final momenta of the proton (P = P' + A). We will be interested in deeply 
virtual Compton scattering where is large compared to the masses and t. Then, 
to leading order in = C • Thus ( plays the role of the Bjorken variable in 

deeply virtual Compton scattering. For a fixed value of —t, the allowed range of ( is 
given by 



< c < 



4M2 , , , 

1 + — - - 1 . (20) 



2M2 \\ (-t) 

The form factor H{x, t, () describes the proton response when the helicity of the proton 
is unchanged, and E{x,t,() is for the case when the proton helicity is flipped. Two 
additional functions H{x, t, (), and E{x, t, () appear, corresponding to the dependence 
of the Compton amplitude on quark helicity. 



Recently, Markus Diehl, Dae Sung Hwang and I have shown how the deeply 
virtual Compton amplitude can be evaluated explicitly in the Fock state representation 
using the matrix elements of the currents and the boost properties of the light-front 
wavefunctions. For the n n diagonal term {An = 0), the arguments of the final-state 
hadron wavefunction are ^jE^, k±i — ^f^A^ for the struck quark and k±i + -^A^ 
for the n — 1 spectators. We thus obtain formulae for the diagonal (parton-number- 
conserving) contribution to the generalized form factors for deeply virtual Compton 



amphtude in the domain [ITS], 0, C < < 1: 
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Figure 2: Light-cone time-ordered contributions to deeply virtual Compton scattering. 
Only the contributions of leading twist in are illustrated. These contributions 
illustrate the factorization property of the leading twist amplitude. 



Vi-C/i )(Xi,t,C) - ^-^===/2(„^„)(Xi,t,C) 



where 

^'i — ) — ^-Li — ^^^± for the struck quark, 

^ . . (23) 

— ■, ^'u — + T^^-i- for ~ 1) spectators. 

A sum over all possible helicities Aj is understood. If quark masses are neglected, the 
currents conserve helicity. We also can check that J2i=i x\ = 1, Z]"=i = 0±. 

For the n + 1 — > n — 1 off-diagonal term (An = —2), consider the case where par- 
tons 1 and n + 1 of the initial wavefunction annihilate into the current leaving n — 1 

— * — * — * 

spectators. Then x^+i — ( — Xi, A;x„+i ~ A± — k±i. The remaining n — 1 partons 

— * 

have total momentum ((1 — C)P'^, — A_l). The final wavefunction then has arguments 

— * — * — * 

— ^^rid /c^j = k±_i -\- ^Aj^. We thus obtain the formulae for the off-diagonal 
matrix element of the Compton amplitude in the domain < xi < C: 



I — 



n+l \ / n+1 



X 



il-n 



{k_Li, k±i, k±n+i = A_L - k±i}, {Ai, A^, A„+i = -Ai}), (24) 



^1 - C ( 1 + 2(^1 Q ) ^ (n+l^n-l) (Xi,t, C) 

^Jo 2(27r)3y 2(27r)3 ^^-^o V 2(27r)3 



n+l \ / n+l 



X 



Hi-E^. 5^^ME^^. [Vi-c]^-" 



'X'4^ln-l)i^iT^'±v ^i) '0(n+l)({^l' ^i' ^"+1 ~ C ~ ^l}) 

{/cxi, k±i, k±n^i — A± — k±i}, {Ai, Aj, A„+i = — Ai}), (25) 
where i = 2,3, ■ ■ ■ ,n label the n — 1 spectator partons which appear in the final-state 
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hadron wavef unction with 



= ^ , ^1. = + Y^^i. ■ (26) 
We can again check that the arguments of the final-state wavefunction satisfy ^27=2 — 

The above representation is the general form for the generalized form factors of the 
deeply virtual Compton amplitude for any composite system. Thus given the light-front 
Fock state wavefunctions of the eigensolutions of the light-front Hamiltonian, we can 
compute the amplitude for virtual Compton scattering including all spin correlations. 
The formulae are accurate to leading order in Radiative corrections to the quark 

Compton amplitude of order asiQ"^) from diagrams in which a hard gluon interacts 
between the two photons have also been neglected. 



5 Applications of QCD Factorization to Hard QCD 
Processes 

Factorization theorems for hard exclusive, semi-exclusive, and diffractive processes al- 
low the separation of soft non-perturbative dynamics of the bound state hadrons from 
the hard dynamics of a perturbatively-calculable quark-gluon scattering amplitude. 
The factorization of inclusive reactions is reviewed in ref. For reviews and bibliography 
of exclusive process calculations in QCD (see Ref. [^, |79t|). 



The light-front formalism provides a physical factorization scheme which conve- 
niently separates and factorizes soft non-perturbative physics from hard perturbative 
dynamics in both exclusive and inclusive reactions pO| . 

In hard inclusive reactions all intermediate states are divided according to A^^ < 
and > domains. The lower mass regime is associated with the quark and gluon 
distributions defined from the absolute squares of the LC wavefunctions in the light cone 
factorization scheme. In the high invariant mass regime, intrinsic transverse momenta 
can be ignored, so that the structure of the process at leading power has the form of 
hard scattering on coUinear quark and gluon constituents, as in the parton model. The 
attachment of gluons from the LC wavefunction to a propagator in a hard subprocess 
is power-law suppressed in LC gauge, so that the minimal quark-gluon particle-number 
subprocesses dominate. It is then straightforward to derive the DGLAP equations from 
the evolution of the distributions with logA^. The anomaly contribution to singlet 
helicity structure function gi{x,Q) can be explicitly identified in the LC factorization 
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scheme as due to the 'y*g qq fusion process. The anomaly contribution would be 
zero if the gluon is on shell. However, if the off-shellness of the state is larger than the 
quark pair mass, one obtains the usual anomaly contribution ||81|| . 

In exclusive amplitudes, the LC wavefunctions are the interpolating amplitudes 
connecting the quark and gluons to the hadronic states. In an exclusive amplitude 
involving a hard scale all intermediate states can be divided according to Ai'^ < 

< and M."^ < A"^ invariant mass domains. The high invariant mass contributions 
to the amplitude has the structure of a hard scattering process Th in which the hadrons 
are replaced by their respective (coUinear) quarks and gluons. In light-cone gauge only 
the minimal Fock states contribute to the leading power-law fall-off of the exclusive 
amplitude. The wavefunctions in the lower invariant mass domain can be integrated 
up to an arbitrary intermediate invariant mass cutoff A. The invariant mass domain 
beyond this cutoff is included in the hard scattering amplitude Th. The Th satisfy 
dimensional counting rules P2[. Final-state and initial state corrections from gluon 



attachments to lines connected to the color-singlet distribution amplitudes cancel at 
leading twist. Explicit examples of perturbative QCD factorization will be discussed 
in more detail in the next section. 

The key non-perturbative input for exclusive processes is thus the gauge and frame 
independent hadron distribution amplitude [^, |[ defined as the integral of the valence 



(lowest particle number) Fock wavefunction; e.g. for the pion 

M^i, A) = 1 d^k^ k^^, A) (27) 

where the global cutoff A is identified with the resolution Q. The distribution amplitude 
controls leading-twist exclusive amplitudes at high momentum transfer, and it can be 
related to the gauge-invariant Bethe-Salpeter wavefunction at equal light-cone time. 
The logarithmic evolution of hadron distribution amplitudes 0i^(xj,Q) can be derived 
from the perturbatively-computable tail of the valence light-front wavefunction in the 
high transverse momentum regime []80| , ^ . The conformal basis for the evolution of the 
three-quark distribution amplitudes for the baryons has recently been obtained 
by V. Braun et al.^. 

The existence of an exact formalism provides a basis for systematic approximations 
and a control over neglected terms. For example, one can analyze exclusive semi- 
leptonic 5-decays which involve hard internal momentum transfer using a perturbative 
QCD formalism ^, |^, |SB|, patterned after the perturbative analysis of form 

factors at large momentum transfer. The hard-scattering analysis again proceeds by 
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writing each hadronic wavefunction as a sum of soft and hard contributions 



^|Jn = ^Pf\Ml < A') + ^^^'\Ml > A'), (28) 

where Ai'^ is the invariant mass of the partons in the n-particle Fock state and A is 
the separation scale. The high internal momentum contributions to the wavefunction 
ipn^'^'^ can be calculated systematically from QCD perturbation theory by iterating the 
gluon exchange kernel. The contributions from high momentum transfer exchange to 
the S-decay amplitude can then be written as a convolution of a hard-scattering quark- 
gluon scattering amplitude Th with the distribution amplitudes 0(xj,A), the valence 
wavefunctions obtained by integrating the constituent momenta up to the separation 
scale Ain < A < Q. Furthermore in processes such as B ^ nD where the pion is 
effectively produced as a rapidly-moving small Fock state with a small color-dipole 
interactions, final state interactions are suppressed by color transparency. This is the 
basis for the perturbative hard-scattering analyses p4|, ^ |8^. In a systematic 



analysis, one can identify the hard PQCD contribution as well as the soft contribution 
from the convolution of the light-front wavefunctions. Furthermore, the hard-scattering 
contribution can be systematically improved. 

Given the solution for the hadronic wavefunctions with Ai'f^ < A^, one can 
construct the wavefunction in the hard regime with A4^^ > A^ using projection operator 
techniques. The construction can be done perturbatively in QCD since only high 
invariant mass, far off-shell matrix elements are involved. One can use this method to 
derive the physical properties of the LC wavefunctions and their matrix elements at high 
invariant mass. Since A4'^ = I]r=i (^ ^'^t"^ ^ ^^^^ method also allows the derivation of 
the asymptotic behavior of light-front wavefunctions at large k±, which in turn leads to 
predictions for the fall-off of form factors and other exclusive matrix elements at large 
momentum transfer, such as the quark counting rules for predicting the nominal power- 
law fall-off of two-body scattering amplitudes at fixed 9cm @ helicity selection 



rules |^8|. The phenomenological successes of these rules can be understood within 



QCD if the coupling av{Q) freezes in a range of relatively small momentum transfer 



6 Two-Photon Processes 

The simplest and perhaps the most elegant illustration of an exclusive reaction in 
QCD is the evaluation of the photon-to-pion transition form factor F^^-^^Q"^) Q 
which is measurable in single-tagged two-photon ee — > eevr^ reactions. The form factor 
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is defined via the invariant amplitude = —ie"^ F.^.y{Q'^)e^'^''^ p'^epQa ■ As in inclusive 
reactions, one must specify a factorization scheme which divides the integration regions 
of the loop integrals into hard and soft momenta, compared to the resolution scale Q. 
At leading twist, the transition form factor then factorizes as a convolution of the 
7*7 —>■ qq amplitude (where the quarks are coUinear with the final state pion) with the 
valence light-front wavefunction of the pion: 

F^m(Q') = ^ £ dxMx, Q)T%Ax, Q')- (29) 

The hard scattering amplitude for 77* qq is T^(x,(5^) = [{1 — x)Q^]~^ x 
(1 + 0{as)) . The leading QCD corrections have been computed by Braaten [0. The 
evaluation of the next-to-leading corrections in the physical ay scheme is given in Ref. 
||89| . For the asymptotic distribution amphtude 0^'^y™P*(x) = A/3/,ra;(l — x) one predicts 
Q^F^AQ^) = 2U (1 - 1^^^^) where Q* = e'^/^Q is the BLM scale for the pion form 
factor. The PQCD predictions have been tested in measurements of 67 en^ by the 
CLEO collaboration 0. See Fig. | (b). The observed flat scaling of the Q^F^^Q^) 
data from = 2 to = 8 GeV^ provides an important confirmation of the applicabil- 
ity of leading twist QCD to this process. The magnitude of Q'^F^^^^Q'^) is remarkably 
consistent with the predicted form, assuming the asymptotic distribution amplitude 
and including the LO QCD radiative correction with av{e~^^'^Q) /n ^ 0.12. One could 
allow for some broadening of the distribution amplitude with a corresponding increase 
in the value of ay at small scales. Radyushkin Ong and Kroll have also 
noted that the scaling and normalization of the photon-to-pion transition form factor 
tends to favor the asymptotic form for the pion distribution amplitude and rules out 



broader distributions such as the two-humped form suggested by QCD sum rules |96 



The two-photon annihilation process 7*7 hadrons, which is measurable in single- 
tagged e~^e~ e"'"e~hadrons events, provides a semi-local probe of C = -|- hadron 
systems 7r°, i]*^, r/c, tt+tt", etc. The 7*7 tt+tt^ hadron pair process is related to 
virtual Compton scattering on a pion target by crossing. The leading twist amplitude 
is sensitive to the 1/x — 1/(1 — x) moment of the two-pion distribution amplitude 
coupled to two valence quarks [|78|, |67| . 



Two-photon reactions, 77 HH at large s = (A;i + ^2)^ and fixed 6*0111, provide a 
particularly important laboratory for testing QCD since these cross-channel "Compton" 
processes are the simplest calculable large-angle exclusive hadronic scattering reactions. 
The helicity structure, and often even the absolute normalization can be rigorously 
computed for each two-photon channel In the case of meson pairs, dimensional 
counting predicts that for large s, s^da / dt^'yy — ^ MM scales at fixed t/s or 6'c.m. up to 
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factors of In s/A^. The angular dependence of the 77 — ^ HH amphtudes can be used to 
determine the shape of the process-independent distribution amphtudes, Q). An 

important feature of the 77 — > MM amphtude for meson pairs is that the contributions 
of LandshofF pitch singularities are power-law suppressed at the Born level - even before 
taking into account Sudakov form factor suppression. There are also no anomalous 
contributions from the x — > 1 endpoint integration region. Thus, as in the calculation 
of the meson form factors, each fixed-angle helicity amplitude can be written to leading 
order in 1/Q in the factorized form [Q^ = Pr = tu/s; Qx = ram{xQ, {I — x)Q)]: 



-^77^MM = _^ dy<pMiy^Qy)THix,y,s,9c.m.^M{x,Qx), (30) 

where Th is the hard-scattering amplitude 77 {(1(i){(1q) for the production of the 
valence quarks coUinear with each meson, and (pMiXjQ) is the amplitude for finding 
the valence q and q with light-front fractions of the meson's momentum, integrated over 
transverse momenta < Q. The contribution of non- valence Fock states are power- law 



suppressed. Furthermore, the helicity-selection rules of perturbative QCD predict 
that vector mesons are produced with opposite helicities to leading order in l/Q and all 
orders in as- The dependence in x and y of several terms in T\^\i is quite similar to that 
appearing in the meson's electromagnetic form factor. Thus much of the dependence 
on (j)M{x, Q) can be eliminated by expressing it in terms of the meson form factor. In 
fact, the ratio of the 77 — tt+tt^ and e"'"e~ amplitudes at large s and fixed 

6cM is nearly insensitive to the running coupling and the shape of the pion distribution 
amphtude: 

f(77--vr+7r-) _ A\F^{s)\' ^^^^ 



f (77-./i+/i-) l-COS^^cm.' 

The comparison of the PQCD prediction for the sum of vr+vr" plus K^K~ channels 
with recent CLEO data | |103|| is shown in Fig. |^. The CLEO data for charged pion and 
kaon pairs show a clear transition to the scaling and angular distribution predicted by 
PQCD for W = \f{s.y-y > 2 GeV. See Fig. |^. It is clearly important to measure 
the magnitude and angular dependence of the two-photon production of neutral pions 
and p'^p~ cross sections in view of the strong sensitivity of these channels to the shape 
of meson distribution amplitudes. QCD also predicts that the production cross section 
for charged p-pairs (with any helicity) is much larger that for that of neutral p pairs, 
particularly at large ^c.m. angles. Similar predictions are possible for other helicity-zero 
mesons. The cross sections for Compton scattering on protons and the crossed reaction 
77 ~^ PP high momentum transfer have also been evaluated [^, providing 
important tests of the proton distribution amplitude. 
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Figure 3: Comparison of the sum of 77 — tt+tt^ and 77 — » K^K~ meson pair produc- 
tion cross sections with the scahng and angular distribution of the perturbative QCD 
prediction pO] . The data are from the CLEO collaboration ||103| . 



It is particularly compelling to see a transition in angular dependence between 
the low energy chiral and PQCD regimes. The success of leading-twist perturbative 
QCD scaling for exclusive processes at presently experimentally accessible momentum 
transfer can be understood if the effective coupling av{Q*) is approximately constant 
at the relatively small scales Q* relevant to the hard scattering amplitudes [Q. The 
evolution of the quark distribution amplitudes In the low-Q* domain at also needs to 
be minimal. Sudakov suppression of the endpoint contributions is also strengthened if 
the coupling is frozen because of the exponentiation of a double logarithmic series. 

A debate has continued |100| , |101| , |102| ] on whether processes such as the pion 
and proton form factors and elastic Compton scattering 7^ — > 7p might be dominated 
by higher-twist mechanisms until very large momentum transfer. If one assumes that 
the light-front wavefunction of the pion has the form V^softla;, = Aexp(— & ^^^^_^^ ), 
then the Feynman endpoint contribution to the overlap integral at small k±_ and x ~ 1 
will dominate the form factor compared to the hard-scattering contribution until very 
large Q^. However, this ansatz for '4'soh{x-,ki_) has no suppression at fcj, = for any 
x; i.e., the wavefunction in the hadron rest frame does not fall-off at all for A;^ = 
and kz —00. Thus such wavefunctions do not represent well soft QCD contributions. 
Endpoint contributions are also suppressed by the QCD Sudakov form factor, reflecting 
the fact that a near-on-shell quark must radiate if it absorbs large momentum. One 
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can show that the leading power dependence of the two-particle light-front Fock 
wavefunction in the endpoint region is 1 — x, giving a meson structure function which 
falls as (1 — x)^ and thus by duality a non-leading contribution to the meson form factor 
F(Q^) oc l/Q^. Thus the dominant contribution to meson form factors comes from the 
hard-scattering regime. 

Radyushkin [p.00|| has argued that the Compton amplitude is dominated by soft end- 
point contributions of the proton wavefunctions where the two photons both interact 
on a quark line carrying nearly all of the proton's momentum. This description appears 
to agree with the Compton data at least at forward angles where — t < 10 GeV^. From 
this viewpoint, the dominance of the factorizable PQCD leading twist contributions 
requires momentum transfers much higher than those currently available. However, 
the endpoint model cannot explain the empirical success of the perturbative QCD 
fixed 9 cm. scaling da / dt{'^p — > TT^n) ~ const at relatively low momentum transfer in 
pion photoproduction ||104|] . 

Clearly much more experimental input on hadron wavefunctions is needed, particu- 
larly from measurements of two-photon exclusive reactions into meson and baryon pairs 
at the high luminosity B factories. For example, the ratio ^(77 — > vr'^vr") /^(77 tt+tt" 



is particularly sensitive to the shape of pion distribution amplitude. Baryon pair pro- 
duction in two-photon reactions at threshold may reveal physics associated with the 
soliton structure of baryons in QCD [p.05| , |106|| . In addition, fixed target experiments 
can provide much more information on fundamental QCD processes such as deeply 
virtual Compton scattering and large angle Compton scattering. 



7 Self- Resolved DifFractive Reactions and Light Cone 
Wavefunctions 

Diffractive multi-jet production in heavy nuclei provides a novel way to measure the 
shape of the LC Fock state wavefunctions and test color transparency. For example, 
consider the reaction [§, |^, ^ nA Jeti + Jet2 + A' at high energy where the nucleus 
A' is left intact in its ground state. The transverse momenta of the jets balance so that 
kj^i + k_i_2 = <?! < A ■ The light-front longitudinal momentum fractions also need to 
add to xi + X2 ~ 1 so that Ap/, < R]^. The process can then occur coherently in the 
nucleus. 

Because of color transparency, the valence wavefunction of the pion with small 
impact separation, will penetrate the nucleus with minimal interactions, diffracting 
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into jet pairs [0. The xi = x, X2 = 1 — x dependence of the di-jet distributions 
will thus reflect the shape of the pion valence light-front wavefunction in x; similarly, 
the k±i — k±2 relative transverse momenta of the jets gives key information on the 
derivative of the underlying shape of the valence pion wavefunction p, |[ |107|| . The 
diffractive nuclear amplitude extrapolated to t = should be linear in nuclear number 
A if color transparency is correct. The integrated diffractive rate should then scale as 
A^/R^ ~ y4^/^. Preliminary results on a diffractive dissociation experiment of this type 
E791 at Fermilab using 500 GeV incident pions on nuclear targets [p.08|| appear to be 
consistent with color transparency. The measured longitudinal momentum distribution 
of the jets [|109|] is consistent with a pion light-cone wavefunction of the pion with the 
shape of the asymptotic distribution amplitude, = -\/3/7rx(l — x). Data 

from CLEO |^ for the 77* — >■ 7r° transition form factor also favor a form for the 
pion distribution amplitude close to the asymptotic solution to the perturbative QCD 
evolution equation [^, |^ . 

The diffractive dissociation of a hadron or nucleus can also occur via the Coulomb 
dissociation of a beam particle on an electron beam {e.g. at HERA or eRHIC) or on 
the strong Coulomb field of a heavy nucleus {e.g. at RHIC or nuclear collisions at 
the LHC) ||107|| . The amplitude for Coulomb exchange at small momentum transfer is 
proportional to the first derivative J2i^ii^'^ of the light-front wavefunction, summed 
over the charged constituents. The Coulomb exchange reactions fall off less fast at high 
transverse momentum compared to pomeron exchange reactions since the light-front 
wavefunction is effective differentiated twice in two-gluon exchange reactions. 

It will also be interesting to study diffractive tri-jet production using proton beams 
pA —>■ Jeti+Jet2+Jet3+y4' to determine the fundamental shape of the 3-quark structure 
of the valence light-front wavefunction of the nucleon at small transverse separation 
0. For example, consider the Coulomb dissociation of a high energy proton at HERA. 
The proton can dissociate into three jets corresponding to the three-quark structure 
of the valence light-front wavefunction. We can demand that the produced hadrons 
all fall outside an opening angle 6 in the proton's fragmentation region. Effectively 
all of the light-front momentum J2j — 1 of the proton's fragments will thus be 
produced outside an "exclusion cone". This then limits the invariant mass of the 
contributing Fock state Ai^ > A^ = P"*"^ sin^ 6'/4 from below, so that perturbative 
QCD counting rules can predict the fall-off in the jet system invariant mass Ai. At 
large invariant mass one expects the three-quark valence Fock state of the proton to 
dominate. The segmentation of the forward detector in azimuthal angle can be 
used to identify structure and correlations associated with the three-quark light-front 
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wavefunction [ |107|| . An interesting possibility is that the distribution amphtude of the 
A(1232) for = 1/2,3/2 is close to the asymptotic form X1X2X3, but that the proton 
distribution amplitude is more complex. This ansatz can also be motivated by assuming 
a quark-diquark structure of the baryon wavef unctions. The differences in shapes of the 
distribution amplitudes could explain why the p ^ A transition form factor appears 
to fall faster at large than the elastic p ^ p and the other p —>■ N* transition form 



factors ||110|| . One can use also measure the dijet structure of real and virtual photons 
beams 'y*A Jeti + Jet2 + A' to measure the shape of the light-front wavefunction for 
transversely-polarized and longitudinally-polarized virtual photons. Such experiments 
will open up a direct window on the amplitude structure of hadrons at short distances. 
The light-front formalism is also applicable to the description of nuclei in terms of 
their nucleonic and mesonic degrees of freedom [ 111 , 112 |. Self-resolving diffractive 
jet reactions in high energy electron-nucleus collisions and hadron-nucleus collisions at 
moderate momentum transfers can thus be used to resolve the light-front wavefunctions 
of nuclei. 



8 Higher Fock States and the Intrinsic Sea 

One can identify two contributions to the heavy quark sea, the "extrinsic" contributions 
which correspond to ordinary gluon splitting, and the "intrinsic" sea which is multi- 
connected via gluons to the valence quarks. The leading l/mg contributions to the 
intrinsic sea of the proton in the heavy quark expansion are proton matrix elements 
of the operator ||113|| rj^Tj^ G a^G pyG""^ which in light-cone gauge rj^A^ = A'^ = 



corresponds to three or four gluon exchange between the heavy-quark loop and the 
proton constituents in the forward virtual Compton amplitude. The intrinsic sea is thus 



sensitive to the hadronic bound-state structure 111 141, |35l|. The maximal contribution of 



the intrinsic heavy quark occurs at Xq ~ m±Q/ J^iJn-^ where m± = \Jm? + A;^; i.e. at 
large xg, since this minimizes the invariant mass A^^. The measurements of the charm 
structure function by the EMC experiment are consistent with intrinsic charm at large 
X in the nucleon with a probability of order 0.6 ± 0.3% [HH] which is consistent with the 



recent estimates based on instanton fluctuations [113 



Chang and Hou [[L15|] have recently discussed the consequences of intrinsic charm 



in heavy quark states such as the -B, A^, and T, such as an anomalous momentum 
distribution for B J/ipX. The characteristic momenta characterizing the B meson 
is most likely higher by a factor of 2 compared to the momentum scale of light mesons. 
This effect is analogous to the higher momentum scale of muonium /i"^e~ versus that 
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of positronium e^e~ in atomic physics because of the larger reduced mass. Thus one 
can expect a higher probabihty for intrinsic charm in heavy hadrons compared to hght 
hadrons. 

One can also distinguish "intrinsic gluons" ||116|| which are associated with multi- 
quark interactions and extrinsic gluon contributions associated with quark substruc- 
ture. One can also use this framework to isolate the physics of the anomaly contribution 



to the Ellis- Jaffe sum rule |^ . Thus neither gluons nor sea quarks are solely generated 
by DGLAP evolution, and one cannot define a resolution scale Qq where the sea or 
gluon degrees of freedom can be neglected. 

It is usually assumed that a heavy quarkonium state such as the J/ ip always decays 
to light hadrons via the annihilation of its heavy quark constituents to gluons. How- 
ever, as Karliner and I ||117|| have shown, the transition J/ip pn can also occur by 
the rearrangement of the cc from the J/ ip into the | qqcc) intrinsic charm Fock state 
of the p or vr. On the other hand, the overlap rearrangement integral in the decay 
ip' piT will be suppressed since the intrinsic charm Fock state radial wavefunction 
of the light hadrons will evidently not have nodes in its radial wavefunction. This 



observation provides a natural explanation of the long-standing puzzle ||118|| why the 
J/ip decays prominently to two-body pseudoscalar- vector final states, breaking hadron 
helicity conservation |]SS|, whereas the ijj' does not. 

The higher Fock state of the proton | uudss) should resemble a | i^A) intermediate 
state, since this minimizes its invariant mass M.. In such a state, the strange quark has 
a higher mean momentum fraction x than the s ||1 19| , H, |12(]|| . Similarly, the helicity of 



the intrinsic strange quark in this configuration will be anti-aligned with the helicity 
of the nucleon ||119| , |120|| . This Q ^ Q asymmetry is a striking feature of the intrinsic 
heavy-quark sea. 



9 Non-Perturbative Solutions of Light-Front Quan- 
tized QCD 

Is there any hope of computing light-front wavefunctions from first principles? The 
solution of the light-front Hamiltonian equation H^^^\ \E') = M^ \ \E') is an eigenvalue 
problem which in principle determines the masses squared of the entire bound and 
continuum spectrum of QCD. If one introduces periodic or anti-periodic boundary 
conditions, the eigenvalue problem is reduced to the diagonalization of a discrete Her- 
mitian matrix representation of H^c^ ■ The light-front momenta satisfy x+ = ^n, and 
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= ^K, where J2i f^i = K. The number of quanta in the contributing Fock states is 
restricted by the choice of harmonic resolution. A cutoff on the invariant mass of the 
Fock states truncates the size of the matrix representation in the transverse momenta. 
This is the essence of the DLCQ method ||121|| , which has now become a standard 
tool for solving both the spectrum and light-front wavefunctions of one-space one-time 
theories - virtually any 1 + 1 quantum field theory, including "reduced QCD" (which 
has both quark and gluonic degrees of freedom) can be completely solved using DLCQ 
122[ |63[| . The method yields not only the bound-state and continuum spectrum, but 



also the light-front wavefunction for each eigensolution ||123| , |124 |. 

In the case of theories in 3+1 dimensions, Hiller, McCartor, and I ||125| , |126|| have 
recently shown that the use of covariant Pauli-Villars regularization with DLCQ allows 
one to obtain the spectrum and light-front wavefunctions of simplified theories, such as 
(3+1) Yukawa theory. Dalley et al. have shown how one can use DLCQ in one space- 
one time, with a transverse lattice to solve mesonic and gluonic states in 3 + 1 QCD 
| 127 |. The spectrum obtained for gluonium states is in remarkable agreement with 
lattice gauge theory results, but with a huge reduction of numerical effort. Hiller and 
I ||128|| have shown how one can use DLCQ to compute the electron magnetic moment 
in QED without resort to perturbation theory. 

One can also formulate DLCQ so that supersymmetry is exactly preserved in the 
discrete approximation, thus combining the power of DLCQ with the beauty of super- 
symmetry [|129| , |130l|131|| . The "SDLCQ" method has been applied to several interesting 
supersymmetric theories, to the analysis of zero modes, vacuum degeneracy, massless 
states, mass gaps, and theories in higher dimensions, and even tests of the Maldacena 
conjecture ||129|| . Broken supersymmetry is interesting in DLCQ, since it may serve as 
a method for regulating non-Abelian theories | p.2(j| |. 

There are also many possibilities for obtaining approximate solutions of light-front 
wavefunctions in QCD. QCD sum rules, lattice gauge theory moments, and QCD in- 
spired models such as the bag model, chiral theories, provide important constraints. 
Guides to the exact behavior of LC wavefunctions in QCD can also be obtained from 
analytic or DLCQ solutions to toy models such as "reduced" QCD{1 + 1). The light- 
front and many-body Schrodinger theory formalisms must match In the nonrelativistic 
limit. 

It would be interesting to see if light-front wavefunctions can incorporate chiral 
constraints such as soliton (Skyrmion) behavior for baryons and other consequences of 
the chiral limit in the soft momentum regime. Solvable theories such as QCD{1 + 1) are 
also useful for understanding such phenomena. It has been shown that the anomaly 
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contribution for the 7r° 77 decay amplitude is satisfied by the hght-front Fock 



formahsm in the hmit where the mass of the pion is hght compared to its size |132 



10 Non-Perturbative Calculations of the Pion Dis- 
tribution Amplitude 

The distribution amphtude (f){x, Q) can be computed from the integral over transverse 
momenta of the renormalized hadron valence wavefunction in the light-cone gauge at 



fixed light-cone time [|T^ 



<P{x,Q)= J d'kle (q^- ^^ j (32) 

where a global cutoff in invariant mass is identified with the resolution Q. The distri- 
bution amplitude 0(x, Q) is boost and gauge invariant and evolves in InQ through an 
evolution equation ^ |^. Since it is formed from the same product of operators 
as the non-singlet structure function, the anomalous dimensions controlling (j){x, Q) 
dependence in the ultraviolet logQ scale are the same as those which appear in the 
DGLAP evolution of structure functions The decay tc fiu normalizes the wave 
function at the origin: ao/6 = dx(j){x^Q) = fT,/{2\/3). One can also compute the 
distribution amplitude from the gauge invariant Bethe-Salpeter wavefunction at equal 
light-cone time. This also allows contact with both QCD sum rules and lattice gauge 
theory; for example, moments of the pion distribution amplitudes have been computed 
in lattice gauge theory |133| , |134 . 



Dalley ||135|| has recently calculated the pion distribution amplitude from QCD 
using a combination of the discretized DLCQ method for the x~ and x~^ light-cone 
coordinates with the transverse lattice method |p.36| , p.37|| in the transverse directions. 



A finite lattice spacing a can be used by choosing the parameters of the effective theory 
in a region of renormalization group stability to respect the required gauge, Poincare, 
chiral, and continuum symmetries. The overall normalization gives f-j^ = 101 MeV 
compared with the experimental value of 93 MeV. Figure | (a) compares the resulting 
DLCQ/transverse lattice pion wavefunction with the best fit to the diffractive di-jet 
data (see the next section) after corrections for hadronization and experimental ac- 
ceptance |jl|. The theoretical curve is somewhat broader than the experimental result. 
However, there are experimental uncertainties from hadronization and theoretical er- 
rors introduced from finite DLCQ resolution, using a nearly massless pion, ambiguities 



32 



in setting the factorization scale Q^, as well as errors in the evolution of the distribution 
amplitude from 1 to 10 GeV^. Instanton models also predict a pion distribution am- 
plitude close to the asymptotic form |138|. In contrast, recent lattice results from Del 



Debbio et al. ||134|| predict a much narrower shape for the pion distribution amplitude 
than the distribution predicted by the transverse lattice. A new result for the proton 
distribution amplitude treating nucleons as chiral solitons has recently been derived by 
Diakonov and Petrov ||139|| . Dyson-Schwinger models [|10| of hadronic Bethe-Salpeter 
wavefunctions can also be used to predict light-cone wavefunctions and hadron distri- 
bution amplitudes by integrating over the relative k~ momentum. There is also the 
possibility of deriving Bethe-Salpeter wavefunctions within light-cone gauge quantized 



QCD [21 1 in order to properly match to the light-cone gauge Fock state decomposition. 
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Figure 4: (a) Preliminary transverse lattice results for the pion distribution amplitude 
at _ 10 GeVl The solid curve is the theoretical prediction from the combined 
DLCQ/transverse lattice method |p.35|| ; the chain line is the experimental result ob- 
tained from jet diffractive dissociation |]l|. Both are normalized to the same area 
for comparison, (b) Scaling of the transition photon to pion transition form factor 
Q'^F^T^o[Q'^). The dotted and solid theoretical curves are the perturbative QCD predic- 
tion at leading and next-to-leading order, respectively, assuming the asymptotic pion 



distribution The data are from the CLEO collaboration 92 



11 Conclusions 

In these lectures I have shown how the universal, process-independent and frame- 
independent light-front Fock-state wavefunctions encode the properties of a hadron 
in terms of its fundamental quark and gluon degrees of freedom. Knowledge of such 
wavefunctions will be critical for progress in understanding exclusive B decays. 
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I have shown how, given the proton's hght-front wavefunctions, one can compute 
not only the moments of the quark and gluon distributions measured in deep inelastic 
lepton-proton scattering, but also the multi-parton correlations which control the dis- 
tribution of particles in the proton fragmentation region and dynamical higher twist 
effects. Light-front wavefunctions also provide a systematic framework for evaluating 
exclusive hadronic matrix elements, including time-like heavy hadron decay amplitudes, 
form factors, and the generalized form factors that appear in deeply virtual Compton 
scattering. The light-front Hamiltonian formalism also provides a physical factorization 
scheme for separating hard and soft contributions in both exclusive and inclusive hard 
processes. 

The leading-twist QCD predictions for exclusive two-photon processes such as the 
photon-to-pion transition form factor and 77 — >■ hadron pairs are based on rigorous 
factorization theorems. The recent data from the CLEO collaboration on F^t^{Q'^) and 
the sum of 77 tt+tt^ and 77 — ^ K^K~ channels are in excellent agreement with the 
QCD predictions. It is particularly compelling to see a transition in angular depen- 
dence between the low energy chiral and PQCD regimes. The success of leading-twist 
perturbative QCD scaling for exclusive processes at presently experimentally accessible 
momentum transfer can be understood if the effective coupling av{Q*) is approximately 
constant at the relatively small scales Q* relevant to the hard scattering amplitudes |]89[ . 
The evolution of the quark distribution amplitudes In the low-Q* domain at also needs 
to be minimal. Sudakov suppression of the endpoint contributions is also strengthened 
if the coupling is frozen because of the exponentiation of a double logarithmic series. 

One of the formidable challenges in QCD is the calculation of non-perturbative 
wavefunctions of hadrons from first principles. The recent calculation of the pion 
distribution amplitude by Dalley [|135|] using light-cone and transverse lattice methods 
is particularly encouraging. The predicted form of (j)T^{x, Q) is somewhat broader than 
but not inconsistent with the asymptotic form favored by the measured normalization 
of Q"^ F^T^o {Q"^) and the pion wavefunction inferred from diffractive di-jet production. 

Clearly much more experimental input on hadron wavefunctions is needed, particu- 
larly from measurements of two-photon exclusive reactions into meson and baryon pairs 
at the high luminosity B factories. For example, the ratio ^(77 — ^ tt'^tt^) /^(77 — ^ tt+tt^) 
is particularly sensitive to the shape of pion distribution amplitude. Baryon pair pro- 
duction in two-photon reactions at threshold may reveal physics associated with the 
soliton structure of baryons in QCD ||105|| . In addition, fixed target experiments can 
provide much more information on fundamental QCD processes such as deeply virtual 
Compton scattering and large angle Compton scattering. 
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A remarkable new type of diffractive jet production reaction, "self-resolving diffrac- 
tive interactions" can provide direct empirical information on the light-front wavefunc- 
tions of hadrons. The recent E791 experiment at Fermilab has not only determined 
the main features of the pion wavefunction, but has also confirmed color transparency, 
a fundamental test of the gauge properties of QCD. Analogous reaction involving nu- 
clear projectiles can resolve the light-front wavefunctions of nuclei in terms of their 
nucleon and mesonic degrees of freedom. It is also possible to measure the light-front 
wavefunctions of atoms through high energy Coulomb dissociation. 

There has been notable progress in computing light-front wavefunctions directly 
from the QCD light-front Hamiltonian, using DLCQ and transverse lattice methods. 
Even without full non-perturbative solutions of QCD, one can envision a program to 
construct the light-front wavefunctions using measured moments constraints from QCD 
sum rules, lattice gauge theory, and data from hard exclusive and inclusive processes. 
One can also be guided by theoretical constraints from perturbation theory which dic- 
tate the asymptotic form of the wavefunctions at large invariant mass, x —>■ 1, and 
high k_i_. One can also use ladder relations which connect Fock states of different parti- 
cle number; perturbatively-motivated numerator spin structures; conformal symmetry, 
guidance from toy models such as "reduced" QCD{1 + 1); and the correspondence to 
Abelian theory for A^^^; — > 0, as well as many-body Schrodinger theory in the nonrela- 
tivistic domain. 
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